Higher Order Corrections to Density and Temperature of Fermions from Quantum 

Fluctuations 
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A novel method to determine the density and temperature of a system based on quantum 
Fermionic fluctuations is generalized to the limit where the reached temperature T is large compared 
to the Fermi energy e/. Quadrupole and particle multiplicity fluctuations relations are derived in 
terms of j- . The relevant Fermi integrals are numerically solved for any values of and compared 
to the analytical approximations. The classical limit is obtained, as expected, in the limit of large 
temperatures and small densities. We propose simple analytical formulas which reproduce the nu- 
merical results, valid for all values of jj. The entropy can also be easily derived from quantum 
fluctuations and give important insight for the behavior of the system near a phase transition. A 
comparison of the quantum entropy to the entropy derived from the ratio of the number of deuterons 
to neutrons gives a very good agreement especially when the density of the system is very low. 

PACS numbers: 25.70.Pq,42.50.Lc, 64.70.Tg 



The availability of heavy-ion accelerators which pro- 
vide colliding nuclei from a few MeV/nucleon to 
GeV/nucleon has fueled a field of research loosely referred 
to as Nuclear Fragmentation. The characteristics of the 
fragments produced depend on the beam energy and the 
target-projectile combinations which can be externally 
controlled [lHl!- Fragmentation experiments could pro- 
vide informations about the nuclear matter properties 
and constrain its EOS Long ago, W. Bauer stressed 
the crucial influence of the Pauli blocking in the mo- 
mentum distributions of nucleons emitted in heavy ion 
collisions near the Fermi energy [5]. We have recently 
proposed a method to estimate the density and temper- 
ature based on fluctuations estimated from an event by 
event determination of fragments arising after the en- 
ergetic collision [6]. A similar approach has also been 
applied to observe experimentally the quenching of fluc- 
tuations in a trapped Fermi gas We go beyond the 
method of 0] by including quadrupole fluctuations as well 
to have a direct measurement of densities and tempera- 
tures for subatomic systems. In this paper, we extend the 
method to derive the entropy of the system and we show 
how to recover the classical limit when the temperatures 
are large compared to the Fermi energy. We apply the 
proposed method to the microscopic CoMD approach [8] 
which includes Fermionic statistics. The resulting energy 
densities and temperatures calculated using protons and 
neutrons display a rapid increase around 3 MeV temper- 
ature which is an indication of a first order phase transi- 
tion. This result is confirmed by the rapid increase of the 
entropy per unit volume in the same temperature region. 
Similar results are found from the entropy density derived 
from the ratio of the number of produced deuterons to 
nucleons. Some differences between the numerical esti- 
mates and the — expansions are found. 
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Figure 1: (Top) — versus multiplicity fluctuations using dif- 
ferent approximations. Full line gives the numerical solution 
of Eq. (j4]), full dots are the lowest order approximation dis- 
cussed in ref.jfjj]; (Bottom) Entropy per particle -S (in units of 
h) versus multiplicity fluctuations. Full line gives the numer- 
ical solution of Eq. iJSJ, full triangles are the Sackur-Tetrod 
results. 



cles. A quadrupole Q xy = {p\ — p^) is defined in a direc- 
tion transverse to the beam axis (z-axis) and the average 
is performed, for a given particle type, over events. Such 
a quantity is zero in the center of mass of the equilibrated 
emitting source. Its variance is given by the simple for- 
mula: 
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A method for measuring the temperature was proposed 
in Q based on momentum fluctuations of detected parti- 



where n(p) is the momentum distribution of particles. In 
Q a classical Maxwell-Boltzmann distribution of parti- 
cles at temperature T c ; was assumed which gives: <j\ v = 
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Figure 2: Fqc versus Symbols as top panel in Fig \T\ 



N4m 2 Tci, m is the mass of the fragment. N is the av- 
erage number of particles. In heavy ion collisions, the 
produced particles do not follow classical statistics thus 
the correct distribution function must be used in Eq. (Q} . 
Protons(p), neutrons(n), tritium etc. follow the Fermi 
statistics while, deuterium, alpha etc., even though they 
are constituted of nucleons, should follow the Bose statis- 
tics. In this work we will concentrate on fermions only 
and in particular p and n which are abundantly produced 
in the collisions thus carrying important informations on 
the densities and temperatures reached. Using a Fermi- 
Dirac distribution n(p) 
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where £ = ^— is the energy , ^ is the chemical potential, 
v = ^. Fqc{v) is the quantum correction which should 
converge to one for high T (classical limit). Expand- 
ing to the lowest order in jj, where Sf = e/o(^) 2 ^ 3 = 
36 (-P )2/3 MeV is the Fermi energy of nuclear matter, the 



following result was obtained in 
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Within the same framework we can calculate the fluctu- 
ations of the p, n multiplicity distributions. These are 
given by [Toj : 
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The lowest order expansion in ( — ), was also derived in 



and is given by: 
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From the above equation (4) we can calculate numer- 
ically the multiplicity fluctuations for a given v and re- 
cover the value of ( ) from the following equation which 
is solved numerically: 
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In Fig. [T] we plot the quantity — vs the normalized 

fluctuations obtained by solving numerically eqs.(4) and 
eqs.(6) while the lowest order approximation, eq.(5), is 
given by the full dots. Since in experiments or modeling 
one recovers the normalized fluctuations, it is better to 
find a relation between the normalized temperatures as 
function of the normalized fluctuations displayed in the 
Fig. 1. It is useful to parametrize the numerical results 
as: 
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which is practically indistinguishable from the numeri- 
cal result (full line) reported in Fig. [1] As expected the 
approximations contained in eq. (5) reproduce the nu- 
merical results (full line) up to rs 0.5. Since from 
experimental data or models it is possible to extract di- 
rectly the normalized fluctuations, one can easily derive 
the value of j- from Eq. ([7]). 

Before proceeding further, it is important to test the 
validity of the approximations for the quadrupole fluctu- 
ations by comparing them to the numerical result solving 
Eq. @ . In Fig. [2] we plot the quantum correction term 
Fqc versus 2L, The difference with the classical case is 

again striking (the Fqc in Eq. ([2|) equal to one for a clas- 
sical perfect gas). For simplicity we can parametrize the 
numerical result with the simple approximation: 
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entropy and fluctuations are strongly correlated [1( 
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Figure 3: Temperature versus density normalized to the 
ground state density po = 0.16/m -3 , derived from quantum 
fluctuations, Eqs. |T][21EJ). Open dots and open squares are 
the approximation at the lowest order in j- , full stars and 

open stars are the classical cases similar to [9J, the full trian- 
gles are the numerical results. The top panel refers to protons 
and the bottom panel refers to neutrons. 



which is indistinguishable from the numerical result dis- 
played in Fig. [2](full line). Clearly such an equation 
converges to one at high T as expected. Notice that the 
lowest order approximation Q is valid up to a modest 
— as 0.2. Eqs.JT]) and KJ| might be very useful when de- 
riving densities and temperatures from data or models, 
without worrying if one is in the classical or fully quan- 
tum limit, the only constraint is that we are dealing with 
fcrmions. 



Once the density and the temperature of the system 
have been determined it is straightforward to derive other 
thermodynamical quantities. One of such quantities is 
the entropy: 
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and z = eT is the fu- 
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U and A are the internal and Helmotz free en- 
ergy respectively [loj . This equation can be numerically 
evaluated and the results are plotted in Fig. QJbottom 
panel). For practical purposes it might be useful to have 
a parametrization of the entropy in terms of the normal- 
ized fluctuations, which is physically transparent since 



S_ 

N 



\fn 



-41.68 



41.68 



/, _ ((AiV)2K ,o2 2 



-2.37 



((AAQ 2 ) 
(N) 



0.83( 



((AiV) 2 ) 
(N) 



f (10) 



The latter fit is indistinguishable from the numerical re- 
sult plotted in Fig. [T] (full line-bottom panel) while the 
Sackur-Tetrod result (full triangles) is valid in the classi- 
cal limit as confirmed in the figure 1. 

To illustrate the strength of our approach we simulated 
40 Ca + 40 Ca heavy ion collisions at fixed impact param- 
eter b = lfm and beam energies E^/A ranging from 4 
MeV/A up to 100 MeV/A. Collisions were followed up to 
a maximum time t = 1000/m/c in order to accumulate 
enough statistics. Particles emitted at later times (evap- 
oration) could affect somehow the results and this might 
be important especially at the lowest beam energies. A 
complete discussion of these simulations can be found in 
@, here we will use the results to compare the different 
approximations. 

In Fig. [3] we plot the temperature vs density as ob- 
tained from the quadrupole and multiplicity fluctuations. 
The top panel refers to protons while the bottom to neu- 
trons. As we can see from the figure, the results obtained 
using the fit functions, Eqs. ([7]) and ((SJ), deviate slightly 
from the lowest order approximations given in Eqs. ^ 
and ([5]). This is a signature that we are in the fully 
quantum regime for the events considered. For compar- 
ison, in the same plot we display the classical temper- 
atures which are systematically higher than the quan- 
tum one, see Eq.([5]) and Fig. [32] Q. We notice that for 
a given excitation energy we can derive a classical or a 
quantum temperature, but the density can be derived for 
the quantum case only within our approach. Of course 
other methods could be devised that give both classical 
temperatures and densities using suitable fragment ratios 
[ID |. We stress that those classical temperatures do not 
need to coincide with the classical temperatures consid- 
ered here since we are dealing with protons and neutrons 
only. Larger fragments could be also included and a dis- 
cussion on this can be found in [9], [l2|, [HI • 

To better summarize the results we plot in Fig. @] (top 
panel) , the energy density s — (^j L )p versus tempera- 
ture Q. Different particle types scale especially at high 
T where Coulomb effects are expected to be small. A 
rapid variation of the energy density is observed around 
T « 2MeV for neutrons and T w 3MeV for protons 
which indicates a first order phase transition [l4| . As we 
see from the figure, the numerical solution of the Fermi 
integrals gives small corrections while keeping the rele- 
vant features obtained in the lowest approximation in- 
tact. This again suggests that in the simulations the 
system is fully quantal. We also notice that Coulomb ef- 
fects become negligible at T >> 3MeV where the phase 
transition occurs. The smaller role of the Coulomb field 
in the phase transition has recently been discussed exper- 



4 



1 

^ 10" 1 

> 

0) 

2 10" 2 

u 

10 3 
10 1 



E 



10 



O A 




O A 



□ T 



O A 



A A 



V 



T(MeV) 



Figure 4: (Top) Energy density versus temperature. Symbols 
as in Fig. (Bottom) entropy density versus temperature. 
The opens symbols refer to the entropy density calculated 
from the ratios of the produced number of deuterons to pro- 
tons (triangles) (neutrons-stars), eq.(ll). 



imentally in the framework of the Landau's description 
of phase transitions (l5| . 

In order to confirm the origin of the phase transition, it 
is useful to derive the entropy density £ = (jj)p which is 
plotted in the bottom panel of Fig. @] The rapid increase 
of the entropy per unit volume is due to the sudden in- 
crease of the number of degrees of freedom (fragments) 
with increasing T. The entropy can be also derived using 
the law of mass action from the ratio of the produced 
number of deuterons to protons (or neutrons) Rd,p(n) 
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The CoMD results from eq.(ll) multiplied by the den- 
sity, are plotted in Figure 4 (bottom panel) with open 
symbols. We find an overall qualitative good agreement 
of the entropy density to the quantum results, eq.10, es- 
pecially for neutrons. Very interesting is the good agree- 
ment for neutrons at low T where the particles are emit- 
ted from the surface of the nuclei which is at low density, 
see also Fig. 3. Such a feature is not present for the pro- 
tons due to larger Coulomb distortions. There is a region 
near the transition (T « 3MeV), where both ratios do 
not reproduce the quantum results. However, at large 
temperatures it seems that all methods converge as ex- 
pected. 

In conclusion, in this work we have addressed a gen- 
eral method for deriving densities and temperatures of 
fcrmions. For high temperatures and small densities 
the classical result is recovered as expected. However, 
we have shown in CoMD calculations that the effect of 
higher order terms give small differences in the physi- 
cal observables considered in this paper but they could 
become large when approaching the classical limit. To 
overcome this problem we have produced suitable param- 
eterizations of quadrupole and multiplicity fluctuations 
which are valid for fermions at all temperatures and den- 
sities. The results obtained in this paper are quite general 
and they could be applied to other systems, for instance 
trapped Fermi gases [7|, to determine the entropy from 
normalized quantum fluctuations. We have also shown 
that the quantum entropy can be compared to the one 
derived from the ratio of the number of deuterons to pro- 
tons or neutrons produced in the collisions. Especially 
the neutrons seem to give cleaner results but of course 
they are more difficult to determine experimentally. 
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